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Abstract. 

Frustrated magnets exhibit unusual critical behaviors: they display scaling 
laws accompanied by nonuniversal critical exponents. This suggests that these 
systems generically undergo very weak first order phase transitions. Moreover, 
the different perturbative approaches used to investigate them are in conflict and 
fail to correctly reproduce their behavior. Using a nonperturbative approach we 
explain the mismatch between the different perturbative approaches and account 
for the nonuniversal scaling observed. 



PACS numbers: 75.10.Hk, 64.60.-i, ll.10.Hi 

1. Introduction 

After twenty five years of investigations, the nature of the phase transition undergone 
by frustrated magnets is still strongly debated ^ 0. The main feature of these 
systems is that, due to frustration, their ground state displays a noncollinear order. 
This is the case for the celebrated Stacked Triangular Antiferromagnets (STA) whose 
ground state exhibits a 120° structure. As a consequence, for N = 3 for instance, 
the £0(3) symmetry group of the high-temperature phase is fully broken in the low- 
temperature one. Also, the order parameter is no longer given by a vector, as in the 
nonfrustrated case, but by a matrix made of two orthonormal vectors 4>\ and (f>2- These 
characteristics have led to a conjecture of the existence of a new universality class for 
frustrated magnets. In fact, the interpretation as well as the description of their critical 
physics appear to be considerably more involved than for nonfrustrated magnets: first, 
because there are several experimental facts incompatible with a universal critical 
behavior; second, because the different theoretical approaches used to analyze their 
behaviors are in conflict and fail to reproduce the experimental data. 

2. Experimental and numerical contexts 

The first indications of a nontrivial critical behavior for frustrated magnets appear 
in the analysis of the experimental and numerical data. An extensive analysis of the 
compounds and systems supposed a priori to belong to the same universality class as 
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the STA, such as Helimagnets, has been performed by the present authors in 02]. 
Without entering into the details we give the main features that have emerged from 
this analysis. First, scaling laws are found in most compounds. This is also the case 
within all — but an important one [S] — numerical simulations of the STA. This 
scaling behavior is associated with critical exponents that differ from those of the 
standard 0(N)/0(N — 1) universality classes. This fact has been at the origin of 
the belief that frustrated magnets could constitute a new universality class. However, 
there are other facts characterizing the critical behavior of both XY and Heisenberg 
systems that go against this belief. 

i) There are two groups of compounds that clearly differ by their critical 
exponents. 

ii) The anomalous dimension 77 is negative for many compounds and for many 
systems numerically simulated. We recall that the anomalous dimension rj is always 
positive if the phase transition is truly of second order and if the underlying field 
theory is a Ginzburg-Landau- Wilson (GLW) </> 4 -like theory, as it is the case here. 

Hi) The scaling relations are violated by several standard deviations for many 
compounds. 

iv) Models differing from STA by their microscopical details undergo either strong 
first order phase transitions (XY case) E] or exhibit scaling behaviors incompatible 
with that of the numerical STA model (Heisenberg case) [3 E] • 

v) Recent Monte Carlo and Monte Carlo RG approaches of STA and GLW models 
give clear indications of first order behaviors [&]• 

From all these observations, it appears that an explanation of the critical physics 
of frustrated magnets assuming a standard second order critical behavior relying on 
the existence of a unique universality class is very likely excluded. 

3. Field theoretical context 

3.1. The nonlinear sigma model approach 

A first theoretical approach consists in a low-temperature expansion performed near 
d = 2 on the suitable field theory, the 0(3) x 0(2)/0(2) nonlinear sigma (NLcr) 
model |l()l 111) which is obtained by expressing the interaction between spins in terms 
of the Goldstone modes. The renormalization of the interaction between these modes 
is treated within a double expansion in the temperature T and in e = d — 2. The 
result of this study is that there exists, above d — 2, a stable fixed point [23 E]- At 
this fixed point, the symmetry breaking scheme is enlarged from 0(3) x 0(2)/0(2) 
to 0(3) x 0(3)/0(3) which itself has the same local structure as (9(4)/0(3) on which 
is built the NLcr model relevant to four-component nonfrustrated spins. As seen from 
the NLcr model approach — that only depends on this local structure — the two 
above theories are equivalent to all orders. Accordingly, the critical properties of 
the STA model between d — 2 and d = 4 should be governed by the standard 
(9(4)/0(3) universality class ^J^J. This is, however, in disagreement with both 
the experimental and numerical data as well as the weak-coupling analysis of the 
GLW model of frustrated magnets around d = 4. Therefore, there must exist 
nonperturbative phenomena — with respect to the expansion parameters T and 
e = d — 2 — leading to the breakdown of the NLcr model approach somewhere between 
d = 2 and d — 3. 
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3.2. The GLW model approach 

An alternative approach is to consider a Ginzburg-Landau- Wilson (GLW) 
representation of the physics of frustrated magnets ^1 El GUI El El- This can 
be obtained by relaxing, via a potential, the constraints defining the ground state. 
The resulting Hamiltonian writes: 




The RG equations for the coupling constants u and v have been first determined in 
a double expansion in u, v and e = 4 — d. At leading order |131 1141 ITS) , they display 
no stable fixed point in the XY and Heisenberg cases. A first order phase transition 
is thus expected in d — 3, in contradiction with the NL<r model results and with the 
existence of scaling behaviors for frustrated magnets. An important generalization of 
the GLW model (QJ is to consider 0i and 02 as A-component vectors ^JEIEI- For 
N larger than a critical value N c (d) depending on the dimension, the RG equations 
display two nontrivial fixed points: one, C+, is stable; the other one, C_, is unstable. 
Above N c (d), the transition is thus predicted to be of second order. As N is lowered 
starting from values of N > N c (d), the two fixed points C+ and C_ get closer and 
finally collapse together for N = N c (d). Below N c (d), there is no longer a stable 
fixed point and the transition is expected to be of first order. The value of N c (d) 
has been computed by several perturbative means: in a double expansion in u, v and 
e = 4 — d up to five-loops |17| and directly in d = 3 in a weak-coupling expansion up 
to six-loops ^HJ- These computations lead to an almost common value N c (d = 3) ~ 6. 
This result could be interpreted as the existence of a first order phase transition in 
d = 3 for both XY and Heisenberg spins. This is only true within the 4 — e expansion. 
In effect — and rather surprisingly — within the six-loop weak-coupling expansion 
performed in d = 3 by Pelissetto et al a fixed point is found for N = 2 and 
N = 3. It has been conjectured by Calabrese et al El that these fixed points 
could explain the spreading of critical exponents encountered in frustrated magnets. 
Indeed, due to the specific structure of the RG flow around the focus fixed point, 
the critical exponents display strong variations along the RG trajectories that could 
explain the lack of universality observed experimentally and numerically. Let us now 
underline several drawbacks of the scenario of Calabrese et al. First, it is based on 
the existence of fixed points that are not related to any already known fixed points. 
In particular, the fixed points found for N = 2 and N = 3 within this computation in 
d = 3 that would govern the critical physics of frustrated magnets are, according to 
Pelissetto et al and Calabrese et al, nonanalytically related to those found in the large- 
N as well as in the 4 — e expansions. This means that there is no way to check their 
existence using these perturbative methods. Second, it is very difficult to account, in 
this framework, for the first order behavior deduced from several numerical simulations 
of XY and Heisenberg systems. Third, it is also very difficult to explain why there 
is no physical system characterized by the exponents associated with the fixed point 
found by Pelissetto et al. Finally, there is no explanation of the breakdown of the 
NLcr model predictions. 
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4. Effective average action approach 

To clarify this situation we have employed a nonperturbative approach based on the 
concept of effective average action [211221, TfcM, which is a coarse grained free energy 
in which only fluctuations with momenta q > k have been integrated out. The field 
4> stands here for an average order parameter at the running scale k. In the limit 
k — > 0, all modes are taken into account and Tk=o identifies with the usual Gibbs free 
energy, i.e. the effective action T. At the scale k = A — the overall cut-off — no 
fluctuation has been integrated and identifies with the microscopic Hamiltonian. 
The ^-dependence of Y k is controlled by an exact evolution equation [231121] : 

, , (2) 

Rk ) 

where t = In k/ A. The trace has to be understood as a momentum integral as well as 
a summation over internal indices. In equation J2J, Tu is the exact field- dependent 
inverse propagator — i.e. the second derivative of Y k — and Rk is the effective infrared 
cut-off which suppresses the propagation of modes with momenta q < k. The effective 
average action I\. is a functional, invariant under the symmetry group of the system 
and thus includes all powers of all invariants — and their derivatives — built from the 
average order parameter. Thus, equation "21 is a nonlinear functional equation, too 
difficult to solve exactly so that 1^ must be truncated. One — fruitful — possibility is 
to perform a double expansion of I\- in derivatives and in fields keeping a finite number 
of monomials in the fields while including the most relevant derivative terms [22 • We 
have first considered the following truncation [3*1 141 125): 




^ (fa • dfa - fa ■ d$X + ^ f ? - Kk) 2 + 



4 V ^ Ti r V 4 V2 V 4 



(3) 



where {uj^i \k, Kk, Pk, Zk} are the coupling constants that parametrize the evolution 
of the model with the scale, while p — Tr l fai> an d t — iTr(*0</>) 2 — ^(Tr t 4>4>) 2 are 
the two independent O(N) x 0(2) invariants built from the average order parameter 
4> = 4>ab, a real N x 2 matrix that gathers the two vectors <f>\ and fa- This ansatz 
realizes the symmetry breaking scheme of frustrated magnets: O(N) x 0(2) — » 
0{N — 2) x 0(2)dia g - We have derived the nonperturbative flow equations for the 
different coupling constants using equations "21 and "JJJ. We have also checked the 
convergence of our results with respect to the order of the truncation of in p and 
r. The explicit recursion equations for the coupling constants arc too lengthy to be 
displayed here ["3EJI2E1 so we concentrate on the physical results. 



4-1. Fixed point structure 

We have numerically analyzed our flow equations by varying the dimension between 
d = 2 and d = 4 for different values of N, identified the different fixed points and 
studied their stability. Around d = 2 we find, for all N > 2, a stable fixed point C + 
which identifies with the NLcr model fixed point. For sufficiently large values of N - 
N > 21.8 — one can follow this fixed point up to d = 4 where the one-loop results of 
the weak-coupling expansion performed around this dimension are recovered. We are 
thus able to smoothly interpolate between the perturbative results obtained in d = 2 



Frustrated magnets in three dimensions: a nonperturbative approach 



5 



and d = 4. On the other hand, for any value of ./V > 2, increasing d, one finds a 
new and unstable fixed point C_ in a dimension di(N). For any value of iV < 21.8 
and as d is further increased, the two fixed points C + and C_ get closer, coalesce and 
finally become complex in a dimension d2(N). Above this dimension d2(N), no real 
fixed point is found, a fact which is interpreted as the sign of a first order behavior. 
This collapse of the fixed points C+ and CL, for different values of N, generates the 
curve N c (d). As already said, for large values of N — typically N of order 21.8 - 
the collapse occurs around d = 4 in agreement with the — standard — weak-coupling 
analysis performed around this dimension. However, for low values of N — N = 3 for 
instance — , this collapse of fixed points operates between d — 2 and d = 3: one has 
d%(N = 3) ~ 2.8. This means that our approach provides a nonperturbative solution 
to the breakdown of the NLc model approach |3l 1251 l2l)j : the stable 0(4)/0(3) fixed 
point obtained within this approach collapses in d ~ 2.8 with another — unstable — 
fixed point which is not seen in a low-temperature expansion. 

4-2. The physics in d = 3 

We now give the results we have obtained in d = 3 within our approach 0j . We 
have computed the value of N c (d) in this dimension and found N c (d = 3) = 5.1. 
Then, we have searched fixed points below this value of N c in order to test the results 
of Pelissetto et al and the conjecture of Calabrese et al that the physics of frustrated 
magnets would be governed by fixed points living in the region N < N c (d — 3). This 
search has been unsuccessful. In fact, within our approach, we have discovered that, 
although there is no fixed point, there exists a whole region in the space of coupling 
constants in which the RG flow is slow. This generalizes the notion of pseudo-fixed 
point introduced by Zumbach |27|. For the systems corresponding to these coupling 
constants, the correlation lengths are very large at the transitions, that are thus very 
weakly of first order. We have found scaling behaviors on large ranges of temperature 
for both XY and Heisenberg magnets (see figure JU). Moreover, we have been able 
to account for the spreading of (pseudo-)critical exponents that characterize these 
systems and, by varying the initial conditions of the RG flow, to reproduce the critical 
exponents of the numerically simulated systems and of many compounds. For XY 
systems, we have found 0.25 < (3 < 0.38 and 0.47 < v < 0.58 and, for Heisenberg 
systems, 0.27 < (3 < 0.42 and 0.56 < v < 0.71, ranges of values in which almost 
all experimental and numerical values enter. Let us give more details. We first 
consider the XY case. According to the experimental values of the critical exponents, 
the compounds split into two groups: STA compounds — CsMnBr3, CsNiCla and 
CsMnI 3 — for which one has, on average, = 0.237(4), i/ = 0.555(21), a = 0.344(5) 
and 7 = 1.075(42), and Helimagnets compounds — Holmium and Dysprosium - 
for which one has (3 = 0.389(7), v = 0.558(25) and 7 = 1.10(5) (a is not given 
since it is poorly determined). Within our computations, we find initial conditions 
of the flow leading to critical exponents close to those of Helimagnets: (3 = 0.38, 
v = 0.58 and 7 = 1.13. We also easily find initial conditions leading to (3 = 0.25, 
corresponding essentially to STA compounds. For such values of (3, we find that v 
varies between 0.47 and 0.49, which is somewhat below the value found for CsMnBr3. 
The same splitting of the compounds happens in the Heisenberg case. There is 
a group composed of Cu(HCOO) 2 2CO(ND 2 ) 2 2D 2 0, Fe[S 2 CN(C 2 H 5 ) 2 ] 2 Cl, VC1 2 and 
VBr 2 characterized by (3 = 0.230(8), a = 0.272(35),^ = 0.62(5) and 7 = 1.105(21). 
There is also a group made of CSN1CI3 , CsMnFj and CsMn(Br 19 ,Io.8i)3 for which 
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Figure 1. Log-log plot of the order parameter m and of the correlation length £ 
for TV = 3 as functions of the reduced temperature t r . 

(3 = 0.287(9), a = 0.243(3), v = 0.585(9) and 7 = 1.181(33). Within our computations 
one can find [3 = 0.27 and v — 0.56 which are close to the experimental values 
of the second group. In contrast, we have not been able to reproduce correctly 
the critical exponents of the first group. This difficulty to reproduce the critical 
exponents originates in the different approximations we have used, in particular the 
truncations in fields and derivatives. Beyond this technical difficulty, the trouble is 
that, once the universality is lost, the pseudo-critical exponents depend on the specific 
compound considered and their precise determination would require to know precisely 
the microscopic structure of the compounds or systems studied — providing the initial 
conditions of the RG flow — and to take into account the full field-dependence of 
the potential in r&. Finally and interestingly, we have also been able to find, in the 
Heisenberg and in the XY cases, initial conditions leading to pseudo-critical exponents 
in good agreement with those obtained within the six-loop calculation JSj. In the 
Heisenberg case, one can find exactly the same values — f3 = 0.30, v = 0.55 and 
7 = 1.06 — than Pelissetto et al I n the XY case, one finds f3 = 0.33, v = 0.56 
and 7 = 1.07 which are very close to the six-loop results: (3 = 0.31(2), v — 0.57(3) 
and 7 = 1.10(4). 

5. Conclusion 

On the basis of their specific symmetry breaking scheme, it has been proposed 
that the critical physics of frustrated systems in d = 3 could be characterized by 
critical exponents associated with a new universality class. From this point of view, 
the study of frustrated magnets has been rather disappointing, the experimental 
and numerical contexts excluding such an hypothesis. At the same time, the 
phenomenology of frustrated magnets has displayed a novel kind of critical behavior: 
the existence of generic scaling without universality 0] . Within the framework of 
our nonperturbative approach, this generic and nonuniversal scaling finds a natural 
explanation in terms of slowness of the flow. This method has also explained the 
mismatch between the different perturbative approaches by means of a mechanism of 
annihilation of fixed points that invalidates the low-temperature perturbative approach 
performed on the NLer model. It would be satisfying to understand the very origin 
of this failure. The influence of nontrivial topological configurations on the phase 
transition in d = 3 has been invoked. It remains however to confirm that these 
configurations indeed play such a crucial role and to understand whether they are really 
responsible for the first order character of the transitions. Finally, and even though 
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it could appear as a rather formal question, it would be also satisfying to clarify the 
discrepancy between the nonperturbative and the six-loop perturbative approaches. 
The matching between the sets of exponents found within these different approaches 
suggests that there exists a common origin to these two sets of exponents. One can 
easily imagine that, while the two kinds of computations would be almost converged 
as for the determination of the exponents themselves, one of these computations would 
not be converged as for the true nature — real or complex — of the fixed points. This 
question will be likely clarified in the near future. 

Note added in proof. Let us mention that a Monte Carlo study j25r has confirmed 
the existence of first order behaviors for a whole family of XY frustrated magnets, in 
agreement with our scenario. 
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